A common magnetic origin for the Invar effects in fee iron-based ferromagnets 
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Using first-principles calculations, in conjunction with Ising magnetism, we undertake a theoretical 
study to elucidate the origin of the experimentally observed Invar effects in disordered fee iron-based 
ferromagnets. First, we show that our theory can account for the Invar effects in iron-nickel alloys, 
the anomalies being driven by the magnetic contributions to the average free energies. Second, 
we present evidence indicating that the relationship between thermal expansion and magnetism is 
essentially the same in all the studied alloys, including those which display the Invar effect and those 
which do not. Hence we propose that magnetism plays a crucial role in determining whether a system 
exhibits normal thermal expansion, the Invar effect, or something else. The crucial determining 
factor is the rate at which the relative orientation of the local magnetic moments of nearest-neighbor 
iron atoms fluctuates as the system is heated. 
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Introduction. Iron-based materials are used in a variety 
of technological applications such as springs in watches, 
car bodies, magnetic sensors, and heads of hard drives. 
Despite their ubiquity in everyday life, they exhibit in- 
triguing phenomena ranging from high-temperature su- 
perconductivity in iron pnictides |1[ and Fermi-liquid 
breakdown in iron-niobium Q to the Invar effect in 
transition- metal alloys [sj. Discovered more than 100 
years ago, Invar iron-based materials display anoma- 
lous thermal expansion over broad temperature ranges. 
Feo.65Nio.35 was among the first in this series of sub- 
stances to be discovered [4]. Subsequently, a number 
of systems were reported, some showing ferromagnetism 
(e.g. Feo.68Pdo.32 [5]) and some antiferromagnetism (e.g. 
FeaTi |). 

Over the years, a paradigm has emerged in which In- 
var behavior in iron-based ferromagnets occurs as a re- 
sult of magnetism. Despite this consensus, the nature of 
the mechanism giving rise to the phenomenon remains 
controversial. A prominent example is the debate over 
whether the anomaly in Feo.65Nio.35 is mainly caused by 
temperature-induced changes in the magnitude of local 
moments [7], temperature-induced changes in the rela- 
tive orientation of neighboring spins [3] , partial chemical 
ordering |8|, non-collinear magnetism [9|, or something 
else [3]. Another open question in the field is: are the 
anomalies observed in other systems such as Feo.72Pto.28 
governed by the same essential physics as that responsi- 
ble for the structural behavior in Feo.65Nio.35 pi, [sHlij ? 

Obviously, a unified theory of thermal expansion in 
iron-based ferromagnets should capture the Invar ef- 
fect in ferromagnetic disordered face-centered cubic (fee) 
Feo.65Nio.35, Feo.72Pto.28, and Feo.68Pdo.32 within a sin- 
gle framework. In principle, the linear thermal expansion 
coefficient, a = {l/a){da/dT)p^ can be derived from the 
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FIG. 1: Linear thermal expansion coefficient of an Fei-xNix 
alloy for a given temperature versus the nickel atomic con- 
centration. Filled diamonds show the results of our method 
at room temperature, and filled triangles at T = 1100 K. The 
corresponding open symbols are experimental data 26]. The 
small filled circles show the results of our method with the 
magnetic contribution to the free energy removed; this un- 
derlines the crucial role of this magnetic contribution for the 
Invar effect. 



configuration- averaged Helmholtz free energy F(a, T). In 
practice, application of density-functional theory (DFT) 
to ab initio calculations of finite-temperature average free 
energies remains difficult, even in the adiabatic approxi- 
mation where electronic, vibrational, and magnetic con- 
tributions are treated separately. Indeed, one of the ma- 
jor issues in implementing this strategy is how to incor- 
porate magnetism correctly. This difficulty may explain 
why only very few previous works [1, [llj which include 
first-principles calculations have explicitly studied ther- 
mal effects on the average length or volume of Fe-Ni, 
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Fe-Pt, or Fe-Pd alloys. 

In this Letter, we investigate theoretically thermal 
expansion and its relation to magnetism in disordered 
fee Fei_a;Niaj alloys with x between 0.35 and 0.8, 
Feo.72Pto.28, and Feo.68Pdo.32. We use an Ising model 
to estimate temperature-dependent magnetic quanti- 
ties such as two-spin correlation functions. On the 
other hand, ah initio total-energy calculations performed 
within the coherent potential approximation (CPA) (l3| . 
and the Debye-Griineisen model ^^], provide comple- 
mentary approaches for determining average free ener- 
gies. 

We believe that our simulations contribute to elucidat- 
ing the essential physics of the Invar phenomenon. A key 
conclusion (with broad implications) is that the Invar ef- 
fects in the Fe-Ni, Fe-Pt, and Fe-Pd alloys may all arise 
as consequences of thermal fluctuations in the relative 
direction of neighboring iron spins. 

Method. To estimate the 'anomalous' contribution to 
the thermal expansion coefficient of an Ye^-x^x alloy as 
a function of temperature, we proceed as follows: 

1. We calculate the average spin-spin correlation func- 
tion of a nearest-neighbor Fe-Fe pair in the system at 
temperature T, {SiSj)FF{T)^ and the analogous quantity 
for A- A pairs, {SiSj)AA{T). Here Si indicates whether 
the Ising spin on site i points up {Si = 1) or down 
{Si = —1). To carry out this point of the procedure, a 
mean-field Ising model of the Miiller- Hesse type jl9| with 
properly chosen exchange constants [20] is employed. 
The ability of an Ising model to deal with magnetism 
in Fe-Ni alloys has been tested in previous work [21|. 

2. We convert the output of our Ising model into two 
quantities: xf^{T) and xa^{T). Here xf^{T) represents 
the concentration of Fe atoms whose spins are up and 
xa^{T) the concentration of A atoms whose spins are 
up, both at temperature T. 

3. We perform ab initio calculations to determine the 
average total energy of an alloy in a ferromagnetic (FM) 
state or a disordered local moment (DLM) state [22l 425j . 
We define this state by fixing the statistics of the local 
moments' orientations to reproduce xf^{T) and XAt(^); 
their magnitudes, however, are allowed to vary. Let us 
call the output of this third step E{xF^{T)^XA^{T)^a)^ 
where a is the average lattice spacing. 

4. We add in a vibrational contribution to the free 
energy, estimated within the Debye-Griineisen model us- 
ing parameters determined from step three. This yields 
a free energy 



F{T,XFt{T),XAt{T),a) = 
E{xFt{T),XAt{T),a) + Fvib(T,a). 



(1) 



5. We minimize the free energy ([T]) with respect to a 
to obtain the equilibrium lattice spacing at temperature 
T, a(T). [a{T) = a(T, x^^(T), XAt(T)) .] 
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FIG. 2: The calculated anomalous contribution (Eq. ([2]), tri- 
angles) and the calculated normal contribution (Eq. Q, cir- 
cles) to the thermal expansion coefficient of an Fei-xNix alloy 
for a given temperature versus the nickel atomic concentra- 
tion. Filled symbols correspond to T = 300 K; open symbols 
correspond to T = HOOK. 



6. We evaluate the anomalous contribution to the ther- 
mal expansion coefficient, which we define as the expan- 
sion that would occur if we changed the magnetic config- 
uration (following the Ising model), but did not otherwise 
heat the system 



aa{T)= lim 



2{r)-a{r,XFr{T),XAr{T)) 
a{T) ST 



where T' = T^ ST. 

The counterpart to the anomalous expansion 
cient aa is the normal expansion coefficient an- 

a{r,XFt{T),XAt{T))-a{T) 



an{T)= lim 



a{T) ST 



(2) 
coeffi- 

(3) 



This normal contribution is the expansion that would 
occur if we heated the system without changing the con- 
figuration of the local moments. 

Results and analysis. Fig. [T] shows a comparison of 
our calculated results with measurements of the linear 
thermal expansion coefficient in Fe-Ni alloys [26|. Our 
relatively simple approach reproduces the experimental 
data strikingly well. The figure also shows the result of 
running the same calculations with the magnetic contri- 
bution to the free energy Fmag removed [27| . In this case 
there is essentially no variation of the thermal expansion 
coefficient with x, which substantiates our claim that the 
magnetic ingredient is crucial to the Invar effect. 

The decomposition of the thermal expansion coefficient 
into its two parts is shown in Fig. O It is clear that 
peculiarities in thermal expansion such as the Invar effect 
arise from the anomalous contribution a^. 

Thus the main interpretive question to ask is: What 
is the magnetic origin of a^? To explore this, consider 
our results presented in Fig. 3 for a number of Fe-based 
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FIG. 3: Positive correlation between the anomalous con- 
tribution to the thermal expansion coefficient (calculated 
using equations ([T]) and ([2|)) and the magnetic quantity 
(1 - xfd{SiSj)FF/dT (calculated using the mean-field Ising 
model), for various Fe-based ferromagnets. 



systems: (i) All of the alloys show a positive correla- 
tion between aa and the product of the concentration 
of nearest-neighbor Fe-Fe pairs with the rate at which 
the average spin-spin correlation function changes with 
temperature [28j , 



(1 



2 d{SiSj) FF 



(4) 



Basically, the larger the rate at which iron moments dis- 
order on heating, the larger the negative deviation of 
thermal expansion from normal behavior, (ii) The slopes 
of the correlation lines are of the same order of magni- 
tude. The existence of a connection between the thermal 
expansion anomaly and the temperature derivative of a 
two-spin correlation function is supported by previous 
analytical work [soj . 

Our analysis indicates that the thermal expansion co- 
efficient is similarly related to the magnetic quantity (4) 
in all the studied ferromagnets (Fei-a^Ni^; with x between 
0.35 and 0.8, Feo.72Pto.28, and Feo.68Pdo.32) [31], and 
that the temperature evolution of the magnetic quan- 
tity (4) plays a major role in determining whether such 
a system exhibits normal expansion, the Invar effect, or 
something else. Features in the structural behavior of the 
systems which were observed experimentally (see the in- 
sets in Figs. 4 and 5), but had remained unexplained, can 
now be interpreted on the basis of the abovementioned 
insight and our theoretical results displayed in Figs. 4 
and 5. (i) The drop in the thermal expansion coeffi- 
cient of Fei_a;Nia; at room temperature when the nickel 
concentration is reduced arises from the steep decrease 
of the magnetic term (1 — x)'^d{SiSj) ff / dT . (ii) The 
fact that the expansivity of Invar Feo.72Pto.28 diminishes 
significantly as T/Tc is raised from 0.2 to 0.5 whereas 
that of Invar Feo.65Nio.35 does not reflects the different 
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FIG. 4: The magnetic quantity (1 — x)'^d{SiSj) ff /dT of an 
Fei-ccNicc alloy at room temperature plotted against the nickel 
concentration, according to the mean-field Ising model. The 
inset shows an experimental result for the thermal expansion 
coefficient of Fe-Ni compounds. 
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FIG. 5: Estimated magnetic quantity of the Invar alloys 
Feo.65Nio.35 (thick solid line), Feo.72Pto.28 (thin solid line), 
and Feo.68Pdo.32 (dashed-dotted line) as a function of the re- 
duced temperature T/Tc, where Tc corresponds to the Curie 
temperature of Feo.68Pdo.32. The inset shows the tempera- 
ture dependence of the thermal expansion coefficient of the 
systems, as determined from experiments 0, [sl]. As a re- 
sult of our mean- field approximation, our theoretical value 
for Tc (782 K) is an overestimate compared to the experimen- 
tal value (637 K) 20]. 



behaviors of their magnetic quantities in this interval: 
(1 — xYd{SiSj) FF / dT decreases drastically in the Fe-Pt 
case, but remains almost constant in that of Fe-Ni. 

To conclude, we point out that it would be highly 
desirable (i) to go beyond the present mean-field Ising 
model; and (ii) to investigate the effect of pressure on 
the temperature-dependence of the magnetic quantities 
[33] . Such investigations may lead to a unified theory for 
the Invar effect and the pressure-induced Invar effect Q 
which has been the subject of recent debate. 
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